SYNCHRONIZATION INDUCED BY INTERMITTENT VERSUS PARTIAL 

DRIVES IN CHAOTIC SYSTEMS 
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We show that the synchronized states of two systems of identical chaotic maps subject to either, 
a common drive that acts with a probability p in time or to the same drive acting on a fraction p of 
the maps, are similar. The synchronization behavior of both systems can be inferred by considering 
the dynamics of a single chaotic map driven with a probability p. The synchronized states for these 
systems are characterized on their common space of parameters. Our results show that the presence 
of a common external drive for all times is not essential for reaching synchronization in a system of 
chaotic oscillators, nor is the simultaneous sharing of the drive by all the elements in the system. 
Rather, a crucial condition for achieving synchronization is the sharing of some minimal, average 
information by the elements in the system over long times. 

PACS numbers: 05.45.-a, 05.45.Xt, 05.45.Ra 



Chaos synchronization has attracted much interest 
from both scientists and engineers by providing insights 
into natural phenomena and motivation for practical ap- 
plications in communications and control [Pecora & Car- 
roll, 1990; Boccaletti et al, 2002; Uchida et ai, 2005; 
Argyris et at, 2005; Pikovsky et at, 2002]. This phe- 
nomenon is commonly observed in unidirectionally cou- 
pled systems, where a distinction can be made between a 
drive or forcing subsystem and another driven or response 
subsystem that possesses chaotic dynamics [Pikovsky et 
ai, 2002]. Complete synchronization occurs when the 
state variables of the driving and the response subsystems 
converge to a single trajectory in phase space. On the 
other hand, generalized synchronization of chaos arises 
when a functional relation different from the identity is 
established between the drive and the response subsys- 
tems [Rulkov et ai, 1995; Abarbanel et ai, 1996; Kapi- 
taniak et ai, 1996; Hunt et ai, 1997; Parlitz & Kokarev, 
1999; Zhou & Roy, 2007]. 

Periodic, chaotic, or stochastic drives have been shown 
to induce generalized synchronization in chaotic systems 
[Maritan & Banavar, 1994; Pikovsky et ai, 2002]. The 
auxiliary system approach [Abarbanel et al., 1996] shows 
that when a response and a replica subsystems are driven 
by the same signal, then the orbits in the phase spaces 
of the response and replica subsystems become identi- 
cal and they can evolve on identical attractors, if their 
initial conditions lie on the same basin of attraction of 
the driven-response system. By extension, an ensemble 
of identical chaotic oscillators can also be synchronized 
by a common drive. The specific functional form of the 
drive is not essential; the basic mechanism that leads to 
synchronization is the sharing of the same information 
by the oscillators for all times. In fact, it is has recently 
been shown that the source of the common influence be- 
ing received by the elements in an extended system is 
irrelevant; it could consist of an external drive, or an 
autonomous global interaction field [Alvarez-Llamoza & 
Cosenza, 20081. At the local level, each element in the 



system is subject to a source that eventually induces com- 
plete or generalized synchronization between the source 
and the element. 

In this paper, we explore another mechanism for syn- 
chronization of a system of driven chaotic elements. We 
consider a system of chaotic elements where the exter- 
nal drive acts intermittently on all the elements with a 
probability p. From the analysis of the dynamics at the 
local level, we extend the auxiliary system approach to a 
situation where the drive is applied only to a fraction p of 
randomly chosen elements in a system. We show that the 
complete synchronized states in both, the intermittently 
and the partially driven systems, are equivalent. Our re- 
sults show that the presence of a common drive for all 
times is not indispensable for reaching synchronization in 
an extended system of chaotic oscillators, nor is the si- 
multaneous sharing of the drive by all the elements in the 
system. Our work is motivated by the practical aspect 
of searching for minimal requirements for the emergence 
of synchronization in dynamical systems. 

We search for minimal conditions for the occurrence 
of synchronization of chaos by using models of coupled 
maps. Let us consider a system of N uniformly, inter- 
mittently driven maps, defined as 



Vi, x'. 



t+i 



s{xl,yt), 



with probability p 

with probability (1 — p) (1) 



Vt+i = 9{.Vt), 



where x\ (i = 1, 2, . . . , N) gives the state of the ith map 
at discrete time t, e is the strength of the coupling to the 
drive g{yt). Each map is subject to the same external in- 
fluence (or lack of it) at any time. The coupling function 
is chosen to have the diffusive form 



s{x\,yt) = {l-e)f{x\) + eg{yt). 



(2) 



We assume a chaotic driven dynamics given by a sin- 
gular map belonging to the family f{xt) = b + \xt\^ , 
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where \z\ < 1, 6 is a real parameter and xt € (— oo,cxd). 
These singular maps exhibit robust chaos, with no pe- 
riodic windows in a finite interval of the parameter b 
[Alvarez-Llamoza et ai, 2008]. Robustness is an impor- 
tant property in applications that require reliable oper- 
ation under chaos in the sense that the chaotic behavior 
cannot be destroyed by small perturbations of the system 
parameters. In particular, we employ the value z = —0.5 
for which the corresponding map displays robust chaotic 
dynamics in the range b S [0.62996, 1.88988]. 

A completely synchronized state in the system Eq. ^ 
is given by xl = xt ~ yt, Vi, and it can occur when g = /. 
On the other hand, if g ^ /, generalized synchroniza- 
tion, characterized by the condition x'l ~ xt yt, may 
also arise in this system for p < 1. A synchronized state 
can be characterized by the asymptotic time-average (a) 
(after discarding a number of transients) of the instanta- 
neous standard deviations at of the distribution of map 



variables xj, defined as 



1 ^ 
N ^ 



1/2 



(3) 



where {xt) is the instantaneous mean of the values x\, Wi. 
Stable synchronization corresponds to (a) = 0. Here we 
use the numerical criterion (cr) < 10^^. Figure[T]shows crt 
as a function of time for the intermittently driven system 
Eq. Il]) subject to different chaotic drives g. 
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FIG. 1: (Jt vs. t for the system of driven maps Eq. ([ij with 
N = 10"* and f{x) = 1.2 + |a:^°■^ for different forms of the 
drive g(yt)- Initial conditions were randomly and uniformly 
distributed such that Xq £ [0, 10]. The drive g{yt) is applied 
with probability p starting at t = 1000. (a) g{yt) = 1.2 + 
|2/t|~°'^, e = 0.7, p = 0.6 (complete synchronization), (b) 
g{yt) = l—2yt, e = 0.8, p — 0.8 (generalized synchronization). 



For a given value of the coupling strength e, there is 
a threshold value of the probability p required to reach 
either type of synchronization. Figure [5] shows the re- 



gions for the complete synchronized states for the sys- 
tem Eq. ^ on the space of parameters {p, e) for different 
orbits of a drive g ~ f ■ In particular, when g = f, 
complete synchronization into an unstable period-m or- 
bit of the map /, defined by /'■'"■' (a;„) = a;„ and satisfying 
n,?=i l/'(^")l > 1; where {xi } are the set of 

consecutive points on this orbit, can also be achieved in 
the system Eq. ([IJ, as shown in Fig.O 
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FIG. 2: Regions for completely synchronized (CS) states 
{{a) < 10^^, after discarding 10* iterates) on the plane (p, e) 
for the system of intermittently driven maps, Eq. ([l| with 
f{x) = 1.2 + \x\-°-^, N = 10*. Dashed line: g{yt) = xi = 
-0.393713; solid line: g{yt) = 1.2 -I- \yt\~°-^; dotted hue: 
g{yt) = {xi = 0.204805, X2 = -1.00968}. The boundaries 
also indicate the onset of stability of those same complete syn- 
chronized states in the partially driven system, Eq. The 
boundaries that separate the stable from the unstable regions 
are given by the corresponding curve = (Eq. (|8}) for the 
driven single map, Eq. Q. 

Figure [3] shows the regions for the generalized synchro- 
nized states of the system Eq. ([T]) on the space of param- 
eters {p, e), with a drive g ^ j ■ 

Since each map in the intermittently driven system 
Eq. ([T]) experiences the same external influence (or none) 
at any time, the properties of this system can be analyzed 
from the behavior of the individual local dynamics. Thus, 
we consider a single, intermittently driven map 



xt^\ 



yt+i 



\ s{.xt,yt), 

9{yt), 



with probability p 

with probabihty {1 — p), (4) 



where xt is the driven or response variable, yt is the drive 
and s(xt, yt) has the same functional form as in Eq. ([2]). 
The auxiliary system approach [Abarbanel et at, 1996] 
implies that a driven map can synchronize on identical 
orbits with another, identically driven map. Thus, the 
occurrence of stable synchronization in the single map 
Eq. (HI) should lead to synchronization in the extended 
system of maps Eq. ([1} , even when the drive acts inter- 
mittently in both cases. 
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The single driven map Eq. ([4]) can be regarded as 
a two-dimensional system. The linear stability condi- 
tion for synchronization requires the knowledge of the 
Lyapunov exponents for such system. These are de- 
fined as Ax = limr^oolnia; and Ay = lim^^oo Inij,, 
where Lx and Ly are the magnitude of the eigenvalues of 

[IltLo^ J(xt, yt)]^^'^ , and J{xt, yt) is the Jacobian matrix 
for the system Eq.(|3]), calculated along an orbit. A given 
orbit {xt, yt} from t~Q to t = T — \ can be separated in 
two subsets, according to the source of the xt variable, 
either coupled or uncoupled to the drive, that we respec- 
tively denote &s A = {{xt.,yt} ■ xt = s{xt^i,yt-i)] pos- 
sessing pT elements, and B = {{xt,j/t} : Xt — /(I't-i)} 
having (1 — p)T elements. We get 



l/T 



.t=o / 



/ n n ^'(^*) ^ \ 

t:XteA t:Xt£B 



l/T 



T-1 



V 



n a'ivt) 



t=0 



(5) 

where Sx = = (1 — e)/'(a;), and K is a poly- 
nomial whose terms contain products of Sx, e, and 
g' {yt) to be evaluated along time. Then Lx = 

\Ylx.eA^^T\x.eBn^t)Y'^ and Ly ^ [YlUa'iyt)]"^- 



Thus we get 



Ar = pin 11— e|+ lim — 



In n i/'(^*)i+in n 

(6) 



(7) 
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FIG. 3: Region for generalized synchronized (GS) states 
{{a) < 10~^, after discarding 10* iterates) for the system 
of maps Eq. ([l}, with f{x) = 1.2 + iV = lO", subject 

to the intermittent drive g{yt) = 1 — 2j/^, on the plane {p, e). 
This same region corresponds to generalized synchronization 
(satisfying A^, < 0) for the single driven map Eq. Q with 
giyt) = l-2yl 



where Ag is the Lyapunov exponent of the map g{yt)- 
Synchronization occurs if the Lyapunov exponent corre- 
sponding to the driven map is negative [Rulkov et at, 
1995]; i.e., A^^ < 0. For a given set of parameter values, 
there is a definite value of the probability p at which the 
exponent Ax changes its sign, from positive to negative, 
signaling the onset of synchronization in the dynamics of 
the two-dimensional system Eq. ([?]). 

Figure [3] shows A^, and Ay as a function of p for the 
driven map Eq. ([4]) for different drives g. 
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FIG. 4: Lyapunov exponents and as functions of p for 
the single driven map Eq. © with f{x) = 1.2+|a;|-° ^ calcu- 
lated over 5 x 10"* iterates after discarding 5 x 10'^ transients 
for each value of p. (a) g{yt) = 1.2-1-12/41""'^, e — 0.6 (complete 
synchronization), (b) g{yt) = 1 — 2j/^, e = 0.8 (generalized 
synchronization) . 

When g = /, the condition Ax < implies complete 
synchronization, where xt ~ yt- In this case we get 



Ay 
Ax 



A/ 

: pln|l 



|l-e| + A/, 



(8) 



where A/ is the Lyapunov exponent of the map /. Fig- 
ure[2]shows the stability boundaries, given by A^ = 0, for 
the completely synchronized states of the system Eq. (|1]) 
on the space of parameters (p, e) for different orbits of a 
drive g{yt). 
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On the other hand, if g f, the condition A^^ < cor- 
responds to generalized synchronization, characterized 
by Xf 7^ yt- Figure [3] shows the stabiUty boundaries, 
given by = 0, for the generahzed synchronized states 
of the system Eq. ([4]) on the space of parameters {p, e) for 
different drives g{yt)- In both, Figs. [2] and [H the curves 
Ax = coincide with the boundaries for stabihty of com- 
plete and generalized synchronization, respectively, in the 
extended system Eq. Thus, the condition Aj, < 
for generalized or complete synchronization in the single 
driven map Eq. (jj]) is equivalent to the condition for the 
stability of the synchronized state = xt, Wi, in the 
intermittently driven system of maps Eq. ([T]). 

The above results are a consequence of the auxiliary 
system approach. Furthermore, the equivalence between 
a single driven map and a system of driven similar maps 
also shows that, under some circumstances, the collective 
behavior of an extended system of interacting elements 
can be inferred by considering the dynamics on a single 
element at the local level. As an application of this idea, 
we consider a partially driven system of maps defined as 



vector form as 



yt+i 



I sixl,yt), 

givt) ■ 



with probability p 

with probability {1 — p) (9) 



The parameter p is the probability of interaction of a map 
with the drive g at a time t. The driven elements are 
randomly chosen with a probability p, so that not all the 
maps in the system receive the same external influence at 
all times. Thus, the average fraction of driven elements in 
the system Eq. ^ at any given time is p. In comparison, 
the forcing of the elements in the intermittently driven 
system Eq. ([T|) is simultaneous and uniform; each map 
receives the same influence from the drive g at any t 
with probability p. 

When system Eq. © gets synchronized, we have — 
Xt- However, the synchronized solution exists only if g ~ 
f. Therefore, only complete synchronization xl = xt = 
yt can occur in this system. 

At the local level, each map in the partially driven 
system Eq. ^ is subject, on the average, to an exter- 
nal forcing g with probability p. Thus, the behavior of 
system Eq. © can also be studied from the behavior of 
the single, intermittently driven map Eq. In partic- 
ular, if the system of maps Eq. ([9]) driven with g = f 
reaches a complete synchronized state for some values of 
parameters, then for this same set of parameters the sin- 
gle driven map Eq. ^ subject to the same drive should 
eventually exhibit a synchronized state similar to that of 
system Eq. Thus, the condition A^, < for com- 

plete synchronization in the single driven map Eq. ([?]), 
that implies stable synchronization in the intermittently 
driven system Eq. ([T]), is also equivalent to the condi- 
tion for the stability of the complete synchronized state 
xl = Xt = yt, yi, in the partially driven system Eq. ([5]). 
To see this, when g = f we express the system Eq. © in 



xt+i = Gff(xt) 



(10) 



where the {N + l)-dimensional vectors xt and f(x() 
have components [xf],; = xl and [f(x4)]i = f{xl), re- 
spectively, for i = 1,...,A^, while [xtj^v+i = yt and 
[f(xO]w+i = giyt). The (iV + 1) X (iV + 1) matrix 
Gf expresses the coupling between the N maps and the 
drive at a time t. The matrix Gj at time t possesses pN 
randomly chosen rows, each having its components Gij 
= 1,2, . . . , {N + 1)) equal to 0, except Gi^i = e and 
Gi_N+i ~ £• The remaining N — pN rows, and the row 
(A^+l), have their components Gij = 0, except Gi,; = 1. 
Thus, at a time t, the matrix G( has the form 



Gt 



/(1-6) ••• 

1 
: '■. 



V 



(1-6) 6 

... 1/ 



(11) 



by 



The long-time evolution of the vector state Xt is given 



xt+i = G f(xt). 



(12) 



where the matrix G is the geometric mean of the prod- 
ucts of all possible configurations of Gt over a long time 
T, 



G 



1/T 



(13) 



We get 



/ (1 - e)P^ 

(1-e)^'^ 



G = 







V 







K\ 
K 



(l-e)P^ K 
1 / 



l/T 



(14) 

where K is a polynomial whose terms contain products 
of (1 — e) and e. 

The linear stability analysis [Waller & Kapral, 1984; 
Kaneko, 1990] of the complete synchronized state /(xj) = 
f{xt) yields 



\ak e 



<1, 



(15) 



where Uk {k = 0, 1 . . . , A^) are the set of eigenvalues of 
the matrix G, with ao = 1 and aj, = (e — 1)^ for fc > 0, 
having iV-fold degeneracy. The eigenvector correspond- 
ing to fc = is homogeneous. Thus perturbations of xj 
along this eigenvector do not destroy the coherence, and 
the stability condition associated with fc = is irrelevant 
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for a synchronized state. The other eigenvectors corre- 
sponding to fc ^ are not homogeneous. Thus, condition 
Eq. (fT5|) with k ^ becomes 



one-dimensional lattice, 



pln|l -e| 4- A/ < 0, 



(16) 



which is the same condition for stability of complete syn- 
chronized states in the single driven map, Eq. 0, when 
g = f ■ Thus, the stability boundary A^,, = in Fig. [2] 
for the driven map with g = f coincides with both, the 
boundary that separates the region where complete syn- 
chronization occurs on the space of parameters (p, e) for 
the partially driven system Eq. ([9]), and the boundary 
for complete synchronization in the intermittently driven 
system Eq. ([1]). However, generalized synchronization 
cannot occur in the former system. 




FIG. 5: Average time Ts to reach complete synchronization 
as a function of the probability p for both, the intermittently 
driven system Eq. ^ (bottom curve) and the partially driven 
system Eq. ([Qj (top curve). For both systems, / — 1.2 — 
|a:|"° '', e = 0.7, = 10**, and g ^ f. The error bars on each 
curve correspond to standard deviations resulting from 100 
realizations of random initial conditions for each value of p. 



Although complete synchronization of chaos in both 
classes of driven systems, Eq. ^ and Eq. ([9|), can be 
characterized from the knowledge of the dynamical re- 
sponse of a single driven map, the transient behavior to 
reach such a state is different in each case. Figure ([5]) 
shows the average time Tg required to attain complete 
synchronization as a function of p in both, a partially 
driven system and an intermittently driven system of 
maps. The times Tg are larger in the first case, i.e., a 
spatially uniform forcing on a system is more efficient 
than a non-uniform one for achieving complete synchro- 
nization. 

Complete synchronization in networks of coupled oscil- 
lators subject to either a uniform drive g with a proba- 
bility p, or to a drive g applied on a random fraction p of 
the elements at all times, may also be equivalent. As an 
illustration, consider a uniformly, intermittently driven 



(1 - f - 7)/(^i) + i [f{xl+') + f{xi-')] + eg{yt), 
with probability p; 

{i-i).f{xi) + -i[fixi+')+f{xr')], 

with probability (1 — p); 

(17) 

where 7 is the local coupling parameter. The analogy can 
be established with a similar lattice subject to a drive 
g{yt) acting on a randomly chosen fraction p of maps, 

(1 - e - j)f{xl) + 1 [f{xl+') + f{xr')] + eg{y,), 
_ with probability p; 

^'+1 = ^ (1 - ^)f(xi) + i [f{xi+') + f(xr')] , 

with probability (1 — p). 

(18) 

Periodic boundary conditions are assumed for both lat- 
tices. Figure [S] shows at as a function of time for both 
systems, Eq. pT)) and Eq. p8)) . for a given example. 
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FIG. 6: at vs. t for (a) the intermittently driven lattice 
Eq. ([T7)l and (b) the partially driven lattice Eq. ^T8\. In 
both cases, f{x) = 1.2 + \x\'°-^ , e = 0.8, 7 = 0.2, p = 0.5, 
N = 10^, and the drive g = f is applied starting at t = 1000. 
Random initial conditions are used. 

In summary, we have shown that the synchronization 
behavior of a system of chaotic maps subject to an exter- 
nal forcing can be inferred from the behavior of a single 
element in the system. The local dynamics can be seen 
as a single driven map; when this drive-response system 
reaches synchronization, the auxiliary system approach 
implies that an ensemble of identical maps subject to 
a similar drive should also synchronize. We have shown 
that the regions of stable complete synchronization in pa- 
rameter space of two systems of identical chaotic maps 
subject to the same drive acting either intermittently in 
time, or partially in space, coincide. At the local level, 
the long-time dynamics of both systems can be character- 
ized by a single chaotic map driven with a probability p. 
In addition, we have shown that complete synchroniza- 
tion in networks of interacting elements subject to either 
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intermittent or partial drives can also be similar. This 
result suggests that these systems possess an ergodicity 
property. 

For the intermittently driven system, the sharing of the 
external influence by the elements takes place only a frac- 
tion of the time. In the case of the partially driven sys- 
tem, the external drive is shared only by a fraction of el- 
ements at any time. Thus, neither the presence of a com- 
mon influence for all times, or the simultaneous sharing 
of the same influence by all the elements, seem essential 
for reaching synchronization in systems of chaotic oscilla- 
tors. What becomes crucial for achieving synchronization 
in both systems is the sharing of some minimal, average 
information by the elements over long times. Future ex- 
tensions of this work should include the investigation of 
some quantity, such as the transfer entropy [Schreiber, 
2000] , for measuring this minimal amount of information 
required for both types of synchronization, and the con- 
sideration of other forms of collective behaviors observed 
in dynamical networks, besides synchronization. 
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